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We consider a system of identical van der Pol oscillators, globally coupled through their 
velocities, and study how the presence of competitive interactions affects its synchronisation 
properties. We will address the question from two points of view. Firstly, we will investi- 
gate the role of competitive interactions on the synchronisation among identical oscillators. 
Then, we will show that the presence of an intermediate fraction of repulsive links results 
in the appearance of macroscopic oscillations at that signal's rhythm, in regions where the 
individual oscillator is unable to synchronise with a weak external signal. 

§1. Introduction 

Synchronisation,® or the ability of coupled oscillators to adjust their rhythms, 
is a property that arises in many systems, from pacemaker cells in the heart firing 
simultaneously as a result of their interaction,® to the fetal heart rate adjusting its 
pace to maternal breathing, as an example of forced synchronisation.® 

Typically, oscillators with different frequencies are able to synchronise due to 
a strong enough positive coupling among units. However, interactions in Nature 
are often repulsive and, surprisingly, it was found that under some particular cir- 
cumstances repulsive interactions can actually enhance synchronisation: thus, the 
presence of negative links can prevent the instability of the fully synchronised state 
when it compensates an excessive number of positive links,® or its sparse presence 
can enhance synchronisation in small- world networks.® Most interestingly - since 
it is not always desirable to achieve a state of full synchronisation - the presence of 
repulsive links can give rise to new forms of synchronisation,® that sometimes can be 
described as glassy or glassy-like.®'®'^®'® Additionally, the beam-forming abilities 
of a system of repulsively coupled Stuart-Landau oscillators were considered inP® 

So far, studies have mostly focused on non-identical phase oscillators, and several 
coupling schemes have been chosen, such as locaP® or long-range,® and purely re- 
pulsive® or assuming a competition between repulsive and attractive?® Like inp® 
we want to isolate the effect of different proportions of repulsive interactions by 
considering identical oscillators. However, rather than establishing how full synchro- 
nisation becomes unstable as the fraction p of repulsive links increases P® our focus 
will be on the characterisation of the different configurations that emerge as p grows, 
and its implications for signal transmission when the system is subjected to an ex- 
ternal forcing. Also, unlike^® we will not consider phase oscillators, but instead van 
der Pol oscillators,^® which implies phase, amplitude and frequency synchronisation 
are taken into consideration. 
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The establishment of the role of the coupling structure on synchronisation, in- 
dependently of the detailed specification of the nodes dynamics, can rely on the 
study of the Laplacian matrixP^'^ We will identify and characterise a transition 
region from synchronisation to desynchronisation by analysing the eigenmodes of the 
Laplacian matrix corresponding to different proportions of repulsive links, adapting 
a formalism developed mP^'® 

The second part of the paper will be devoted to explore the role of competitive 
interactions in the synchronisation of the system with an external periodic signal. 
We will choose a signal whose frequency lies outside the region of entrainment for an 
uncoupled oscillator, as well as for an all attractively coupled system. This problem 
is closely related to a second theoretical framework, that of resonance studies, that 
emphasise the importance of an intermediate disorder on the response to an weak 
signal, where disorder can be noisepD'EZJ diversitypS'^'^ or competitive interac- 
tions.^'^ 1 In the latter cases^'® it was found that an intermediate fraction of 
repulsive links was able to amplify the response to an external signal, in bistable sys- 
tems where the external signal was be the only source of movement. In the present 
case, an optimal response should correspond to an adjustment between the intrin- 
sic frequency of the oscillators and that of the external signal; as we will see, that 
optimal response is achieved at an intermediate proportion of repulsive links in the 
case of strong fast signals, whereas weak slow signals are best responded when all 
the links are negative. 

The outline of this paper is as follows: in section [2] we will introduce the model; 
we show that an increase of the proportion of repulsive links leads to loss of synchro- 
nisation in section [3l and in section [J] we show how the presence of repulsive links 
accounts for an enhanced response to external signals; in section [5] we will briefly 
mention some extensions; and conclusions are drawn in section [6j 



§2. Model 

We consider an ensemble of van der Pol oscillators^ { Xi (t), i = 1, . . . , N}, glob- 
ally coupled through their velocities Xi, and subjected to an external periodic forcing 
of amplitude A and frequency fi. The dynamics is described by: 

C N 

x\ = -Xi + - xf)xi + — ^2 Jij (xj - Xi) + Asm(f2t), (2-1) 

where the nonlinearity parameter fi is a positive constant and C is the coupling 
strength. 

The coupling between the oscillators i and j is given by the interaction term 
and can be attractive or repulsive, according to a given probability p. 

j = j = f-1, with probability p, 

1 1, with probability 1 — p. 

The single van der Pol oscillator is a paradigmatic example of a non-linear os- 
cillator. It possesses a stable limit cycle as a result of its nonlinear damping term 
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p(l — xf): for small oscillations, < 1, the system experiences negative damp- 
ing and the oscillations grow, while for \x%\ > 1, the positive damping causes the 
oscillations to decrease. Therefore, independently of the initial conditions or small 
perturbations, its amplitude of oscillations reaches a constant (equal to 2), while 
its detailed shape and period T depend on p, approaching T « (3 — 21n2)/x for 
large p. In this case of large \i S> 1, the oscillations are called relaxational and are 
characterised by the presence of discontinuous jumps intercalated by periods of slow 
motion. 

§3. Desynchronisation among the unforced oscillators, A = 

Fig-fflshows the trajectory (left panels) and respective limit cycles (right panels) 
of two typical individual oscillators for some probabilities p of repulsive links. In 
all cases, the essential characteristics of the van der Pol oscillator, such as a steady 
amplitude and the existence of two time scales, are preserved by this type of coupling, 
as it is reflected by the fact that the stable limit cycles (Fig. [H right panels) maintain 
its basic shape. Having identical natural frequencies, when the coupling constant C 
is strong enough the position of the oscillators becomes synchronised when all the 
interactions are attractive. (Fig. [IJ for p = 0). As the proportion of repulsive links 
grows, both the amplitude and the phase of the oscillators start to desynchronise 
(Fig. [H p = 0.4). Finally, a further increase in the proportion of repulsive links 
(Fig. [U p = 0.60 and p = 1.0), drives the system to a configuration where the global 
variable X(t) = jjYli x i(t) ^ s zero, with several groups oscillating in anti-phase, 
with a decreased frequency and an increased amplitude of oscillations: if we consider 
two groups with the same size that are oscillating in anti-phase, it can be shown 
that their amplitude becomes rescaled as a result of their interaction by a factor of 

We can describe the last configuration characterised by a zero value of the av- 
erage position X(t) as a disordered situation at the macroscopic level. To quantify 
this disordering role of repulsive links, we define, followingpSJ the complex variable 
— I ^'•^it th.6 average z — X^i=i %i the variance of z% normalised by the 
average value of the modulus squared, a 2 [ Zi ]: 

a[Zi] = \ tf-EjLxN' /' (3 " 1} 

here and henceforth (• • • ) denotes a time average. 

The normalised variance can take values between = 1 for maximum disor- 

der, and cr 2 [zi] = when all oscillators are synchronised amongst themselves. From 
this, we choose^ 1 a measure of order that reduces to the Kuramoto order parame- 
ter^ w h en a r| un itg oscillate with the same amplitude: 

p = aA-^M- (3-2) 
As dispersion increases, p decreases from p = 1 to p = 0. As we show in Fig 
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Fig. 1. Trajectories (left panels) and phase portraits (right panels) of two individual oscillators, 
for various probabilities p of repulsive links. Regardless of initial conditions, the van der Pol 
oscillator reaches steady state oscillations with a constant amplitude. In the right panel we 
represent the corresponding limit cycles in the phase plane of (x, x): as p grows the distinction 
between slow and fast motion becomes clearer which is manifested in the more abrupt angles in 
the limit cycles for p = 0.6 and p = 1.0. N=100, C=20, ^ = 10, A=0. 



[21 there is a clear synchronisation-desynchronisation transition for an intermediate 
fraction of repulsive links, that does not depend much on the coupling strength C. 

To characterise desynchronisation further, it is useful to look into the behaviour 
of the field that an oscillator feels as a result of the interaction with other units. 
The average number of effective links F = Yltji Jij = j^i^P ~ 1)C^ ~~ 1)j m a 
given run can in general be different from the particular number an oscillator has, 
fi = ji^2iL±Jij, given that the probability of repulsive links p follows a binomial 
distribution with the corresponding variance. We want to know if there is a correla- 
tion between the fraction of repulsive links an oscillator has and its synchronisation 
with the overall majority. That is described by the following quantity G: 

G = \}j X>v^ [/* - F y = y*I> ifi - F \J > ( 3 - 3 ) 

being X = jj Xj the mean velocity. We observe that the order-disorder transition 
region p ~ [0.4, 0.45] we identify in the left panel of Fig. [21 is accompanied by an 
increase in the influence on an oscillator of its particular coupling configuration, as 
signalled by the peak in G (Fig. [21 right panel). The oscillators with a higher than 
average number of repulsive links form a loosely synchronised group in a different 
slow region than the one where the majority concentrates. 

The partial independence of the state of the oscillator on the global configuration 
opens the possibility of the existence of several different global states, and thus 
hints at the existence of metastable states. We can relate this behaviour to the 
coupling structure, by the spectral analysis of the associated Laplacian matrix J-j = 
Jij — o~ij ^2k=i JkjW^ where 5ij is the Kronecker's delta. We begin by rewriting Eq. 



Synchronisation Induced by Repulsive Interactions 



5 




P P 

Fig. 2. The measures of disorder p (left panel) and G (right panel). Averages over 100 runs, and 
N, C, n and A as in Fig. [1] 



(|2-ip as system of two equations that highlight a fast motion for the Xi variable and 
a slow motion for the yi variable: 



x i y< 



Xi = n 
Vi = - [xi - Asm(Qt)\ 



n N 



ij X 3 



3=1 



(3-4a) 
(3-4b) 



c 



where D 

We focus on Eq. (|3-4ap . letting the slow variable yi be a constant. The Laplacian 
appears naturally in the coupling term, and we can see its positive eigenvalues should 
help deviations from a given state along the x variable We now introduce the 
eigenvalues Q a and eigenvectors e a = (e", . . . , e^r) of the Laplacian matrix, with the 
normalisation condition Y2i e ? e i = 



N 

E 



j'. P a 

l 3 3 



Qa&i 



(3-5) 



Let's us assume the state of a unit i is x° at a given time, where x° is drawn 
from any symmetric random distribution. We perturb the initial states as x° + Sj, 
and express Sj in terms of the eigenvalues and eigenvectors of the Laplacian, so that 
Si = Yla=i Batf- We aim to see how the interaction with other units influences the 
reaction to perturbations. 

We will assume for simplification that x°ef ~ S»=i x i Sj=i e ? = 0> an d 

that the modes are uncoupled. After expanding the equation in terms of the Lapla- 
cian eigenvalues and eigenvectors, we then multiply the resulting equation by ef and 
average over all elements, to find the evolution for the amplitude of the a-th mode: 



G 
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dB, 



dt 



a 



3^ + PR, 



-a 



( 




(3-6) 



where k = jj ~Yld = \{x1) 2 — 1 is a quantity related to the variance of the initial 
conditions, and the participation ratio PR a = l/X]^Li[ e ?] 4 ^ s a classical measure 
of localisation that estimates the number of oscillators that participate significantly 
in a state e a : for a state localised on a fraction / of elements, PR Q tends to /. 
According to Eq. (|3-6p . unless Q a > the amplitude of the mode B a tends to 
zero, and any deviation from the initial state vanishes. Otherwise, mode a is said to 
be an open mode. 

In a precise way, we will define "localised" modes as the ones whose participation 
ratio is less than 0.1 AT, and define a measure M of localisation^ as M = j^jj , where 
Nl is the number of open localised modes, i.e. those satisfying PRq, < O.liV, and 
No is the total number of open modes a. While for extreme probabilities p of 
repulsive links the number of possible values for Q a is very restricted, we find that 
for intermediate levels of p the distribution of Q a is broader and the eigenvalues at 
both tails of the spectrum are localised: the left panel of Fig. [3] illustrates this fact 
for some examples of p. 

In our regime of strong non-linearity, or fj, S> 1, the time an oscillator spends 
on a fast motion is very close to zero; in such case we should have k ~ and the 
condition for open modes can be fulfilled at low enough probabilities of repulsive 
links. Of the open modes, those that correspond to positive Laplacian eigenvalues 
facilitate the growth of perturbations, while those that are negative inhibit it. We 
wish to identify a probability of repulsive links that makes the system flexible enough 
to sustain deviations from the initial symmetric configuration, yet stable enough so 
that perturbations don't spread immediately throughout the entire system, as that is 
what we observe at the transition region (Fig. [2] right panel), where some oscillators 
are more loosely synchronised than others. We can anticipate^ that such situation 
corresponds to a localisation^'^} of positive Laplacian modes. 

As we see in Fig. [3l the peak in G, that signals the range of p where there is an 
heightened dependence of the state of an oscillator on its coupling structure (Fig. 
[2]), coincides with a localisation of the positive eigenvalues of the Laplacian. Those 
localised positive open modes are responsible for keeping the loss of synchronisation 
within controllable limits. 



In this section, we will see how competitive interactions affect the response 
to an external periodic signal. Since in general there can be several frequencies 
present in the output of the global variable X(t) = j? ^ Xi(t), we say the system is 
synchronised with the external signal when the highest peak in the Fourier spectrum 
corresponds to that frequency. 

When the natural frequency of oscillations coincides with the external forcing 
frequency, synchronisation is achieved for vanishing A, and as the two frequencies 



§4. Synchronisation with the external signal, A^O 
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Fig. 3. Left panel: we plot the participation ratio PR a for the Q a eigenvalues. As the probability of 
repulsive links p grows the eigenvalues become dislocated towards higher values. The eigenvalues 
whose participation ratio is below the dashed line are localised. Right panel: the measures of 
localisation M, the measure of disorder G, and the spectral power amplification R, all have the 
maximum value at roughly the same p. For better viewing, we multiplied M by 4, and G by 0.1 
N=100. In the case of M and k=0.2 we have C=20 and fj, = 10. 

diverge, stronger forcings are needed to entrain the system. We will call a signal 
strong when its amplitude is greater than the amplitude of oscillation of the unforced 
van der Pol unit, and we will call it fast when its frequency is higher than the natural 
frequency of the individual van der Pol oscillator. 

We will distinguish between strong fast and weak slow signals, because the mech- 
anism of synchronisation differs in the two cases, although, in both cases, competitive 
interactions are required for an enhanced response. 

4.1. Strong fast signals benefit from intermediate p 

In Fig. U] we plot the synchronisation regions and their relative strength, as 
measured by the spectral power amplification factor^ _R, given by: 



R is roughly proportional to the square of the normalised amplitude of the oscillations 
of X(t) at the frequency Q, being R < 1 when the amplitude of oscillations of the 
forced system is smaller than the amplitude of the external signal. 

When p = 0, (Fig. 31 left panel) the synchronisation region with respect to 
the frequency Q and amplitude A of the external signal has the typical triangular- 
like shape seen on Arnold tonguesP An intermediate fraction of repulsive links 
(p = 0.43, Fig. [H right panel) pushes the synchronisation borders beyond the p = 
values, allowing for synchronisation of faster signals at weaker forcing. 

Fig. [5] shows the steady-state trajectory of the macroscopic variable X(t), for 
different probabilities p of repulsive links, and illustrates the fact a certain proportion 
of repulsive links is required for the system to adjust its rhythm to that of the external 
signal (Fig. p = 0.40), whereas Fig. [6] confirms that this optimal response only 
occurs for an intermediate range of the probability of repulsive links. It should be 
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Fig. 4. We plot the spectral power amplification R in the synchronisation regions for p — (left 
panel) and p = 0.43 (right panel). For better viewing, we use a color codes that saturates for 
R > 2. N, C and /j, as in Fig. Q] 



noted that when entrained, the oscillators adjust their frequency while keeping their 
natural amplitude (Fig. [5]), therefore, the spectral amplification factor R is smaller 
than 1, when A > 2. As expected, the more the natural frequency deviates from the 
forcing frequency, the stronger the signal needs to be in order to entrain the system: 
namely (Fig. [6]), for a forcing frequency Q = 1, the signal strength needs to be 
A = 12 instead of A = 5, when the natural frequency to = 2n/T is ~ 0.19 (/U = 20) 
instead of 0.39 (fi = 10). 
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Fig. 5. Time evolution of the macroscopic variable X(t) when the system is forced by an external 
sinusoidal fast signal (lighter color) of amplitude A = 5 and frequency Q = 1.0, for several 
probabilities of repulsive links p. N, C and /x as in Fig. \T\ 
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Fig. 6. Spectral power amplification, for C = 20, Q = 1 and several system sizes N, averages over 
100 runs. We note that smaller systems become synchronised at lower fractions of repulsive 
links, and are not so dependent on the precise fraction of repulsive links. Additionally, we also 
observed (figures not shown) a resonance with system size for different probabilities of repulsive 
links. This kind of dependence has been explained elsewhere.^ 



To understand the significance of competitive interactions we recall the results 
of the last section. The probability region where weak fast signals can be entrained 
is signalled by the peak in the spectral power amplification R, and coincides with 
the localisation region, as given by the peak in M (Fig. [3]). 

This localisation, or controlled disorder, is crucial for an enhanced response to 
fast signals, for it allows the amplitude of oscillations to be deformed, varying slightly 
enough to place some oscillators closer to the fast motion region, thus allowing a jump 
to another slow region under an external forcing; but not so much that there is a risk 
they would trigger a chain reaction, unless there is a permanent pulling into some 
direction - that is, unless the signal is acting. 

4.2. Weak slow signals benefit from very high p 

In the previous section, we chose to measure the enhancement at the collective 
level, using the macroscopic variable X(t) in our measure of response R, Eq. (|4-1|) ; 
that corresponded to a synchronisation with the external forcing at the individual 
level: the greater the number of entrained oscillators, the greater the response was. 

We find a different situation when we subject our system to a very weak slow 
signal with an amplitude that is smaller than 1, say A = 0.9. A complete ampli- 
tude and frequency synchronisation with this forcing would imply a fast motion, in 
the interval [—1, 1], without any intercalating period of slow motion, thus basically 
destroying the defining feature of a relaxational oscillator (Section [2]). We find it 
impossible for such a weak signal to entrain an individual oscillator. And yet, we 
observe that for a high enough fraction of repulsive links (insets Fig. [7] and Fig. [HJ), 
there is a nearly coincidence between the trajectory of the global variable X(t) and 
the forcing As'm(f2t), with an almost imperceptible phase delay. Therefore, the sim- 



10 



Teresa Vaz Martins, Raul Toral 



plest measure of entrainment, that falls to zero if there is a perfect synchronisation, 
is: 



D = ~ T^77yT\ > (4-2) 



[X(t) -^sin(^)]^ 

The results plotted in Fig. [7] show how the synchronisation with the external signal 
as the fraction p of repulsive links mirrors the loss of synchronisation seen in Fig. [2] 
for the unforced system. 

Again, the mechanism has its roots on the disorder induced by the presence 
of repulsive links. As we saw, after the transition region (Fig. [21 p ~ 0.45), the 
unforced system tends to adopt a configuration corresponding to a zero value of the 
macroscopic variable X (Fig. [5j p = 0.6 and p = 1.0). On the other hand, the forcing 
induces an asymmetry in the oscillations favouring the time spent on the side of the 
signal's extremum, as seen in the right panel of Fig. [HJ Slight as this asymmetry 
may be, it is enough to cause the superposition of the individual waves on the same 
side as the signal's extremum, while on the other side, the oscillations cancel each 
other (see right panel of Fig. [8] for p = 1.0). The total synchronisation of the system 
in case p is zero or small, naturally prevents this phenomenon to happen (see right 
panel of Fig. [8] for p = 0.0). 




§5. Further applications: FitzHugh-Nagumo 



The single uncoupled van der Pol oscillator can be transformed either into a 
linear oscillator by taking fi = or by replacing the nonlinear damping term //(l— x?), 
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Fig. 8. The slower oscillation corresponds to the external signal while the higher frequency os- 
cillations correspond to either the trajectory of the macroscopic variable X(t) or two typical 
individual trajectories. Upper panels: When p = 1.0 we observe the synchronisation of the 
macroscopic variable with a signal that is very weak and slow. Lower panels: we zoom, and plot 
two representative individual trajectories. A = 0.9, Q = 0.01. N, C and fi as in Fig. [T] 



Eq. (|2-1|) by a constant; or into an excitable system - a simplified FitzHugh-Nagumo 
- by adding a constant a such that \a\ > 1 to Eq. (|3-4bp . So, a first direct extension 
consists in a brief exploration of how those transformations affect our results. 

Not surprisingly, we didn't find an enhanced response for linear oscillators: both 
mechanisms of enhancement for slow and fast signals rely on the existence of a 
slow motion region. This situation contrasts with the case studied inP3 In that 
paper, the authors studied a system of linear oscillators with a distribution of natural 
frequencies. Defining as a measure of diversity the variance of the natural frequencies, 
they found an optimal response to an external signal for an intermediate level of 
diversity. Interestingly enough, the enhancement of response also had its origins in 
an intermediate level of disorder. However, the microscopic mechanism was rather 
different: some oscillators had a natural frequency that resonated with the signal's 
frequency, and where able to pull the others due to the positive coupling. In our 
case, there isn't any single oscillator whose frequency can be entrained by the external 
signal. 

On the other hand, the mechanisms we proposed should be applicable to the 
FitzHugh-Nagumo model. Adding a constant a to Eq. (|3-4b|) the system becomes: 



3 Xi 



D N 

y* + n E 4 



ij x j 



(54a) 



i=i 
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Vi = - [xi -Asm(f2t) + a] . (5-lb) 



The interaction via competitive interactions can play the same role of noise or 
diversity, thus enabling rhythmic excursions away from the fixed point. The result 
shown in Fig. [9]bears some similarity with the phenomenon by which we observe that 
the periodicity of oscillations becomes maximally ordered for an intermediate level of 
noise EDJM^ESD diversity,^ or competitive interactions!^ 1 '^ 1 In our case, however, 
and like it was observed inP3 for the case of active rotators, we don't observe any 
oscillations at all unless some interactions are repulsive. 

When we force the excitable system by a strong enough fast signal (Fig. [10]) . 
it starts to oscillate even for p = 0, and for an intermediate amount of repulsive 
interactions the main frequency of oscillations coincides with the external signal 
(Fig. Mp = 0.40). 

On the other hand, when we force a system of FitzHugh-Nagumo elements by 
slow weak signals (Fig. [TTj) we observe bursts with the periodicity of the signal 
for p = 0, while for p = 1 the global variable roughly oscillates along with the 
external forcing. Even though the periodicity of the external signal is detected for 
all fractions of repulsive links, we can imagine situations where we actually want to 
replicate the behaviour of the external signal, and that is only possible when the 
fraction of repulsive links is large enough. 

Both of these results are expected taking into account the arguments we gave 
for the van der Pol oscillator case. 
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Fig. 9. The trajectory of the global variable X for different fractions of repulsive links p in the un- 
forced FitzHugh-Nagumo system showing a similar phenomenon to coherence resonance. Other 
parameters: a = 1.1, N, C and fx as in Fig. [T] 
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Fig. 10. Trajectories of the global variable X for a system of FitzHugh-Nagumo units in the ex- 
citable regime forced by a fast signal, for increasing fractions of repulsive links p. Other param- 
eters: a = 1.1, A=12, Q = 1, N, C and as in Fig. [1] The pointed line shows the external 
signal multiplied by 0.1, for better viewing. 
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Fig. 11. Representative trajectories of the global variable for a system of FitzHugh-Nagumo units 
in the excitable regime, a — 1.1, for some probabilities p of repulsive links. Other parameters: 
Q — 0.01, A=0.9, and N, C and fj, as in Fig. [T] The signal is represented in a lighter color. 



§6. Conclusions 

We have shown that the presence of repulsive links in a system of globally cou- 
pled van der Pol oscillators can enhance the response to an external signal. This 
phenomenon is verified regardless of whether the signal is strong and fast, or weak 
and slow, and it is in every case directly related to a loss of synchronisation and the 



14 



Teresa Vaz Martins, Raul Toral 



existence of a slow motion region, but the microscopic mechanism of enhancement 
is different in each case. 

From the point of view of a strong fast signal, the van der Pol oscillator can be 
approximated by a bistable system, implying a threshold that is regularly overcome 
with the help of an intermediate proportion of repulsive links, by means of the de- 
formation of the slow motion region. In the case of very slow signals, the mechanism 
is associated with the tendentiously anti-phase oscillations that occur when there is 
a majority of repulsive links. 

In both cases, enhancement is directly related to a loss of full synchronisation 
when the fraction of repulsive links increases. One can imagine that a different 
coupling scheme might enhance or hinder the results we found, since it is known 
that the network topology plays a role in synchronisation properties.^ 
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